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We examine the molecular-atomic transition in liquid hydrogen as it relates to met- 
allization. Pair potentials are obtained from first principles molecular dynamics and 
compared with potentials derived from quadratic response. The results provide in- 
sight into the nature of covalent bonding under extreme conditions. Based on this 
analysis, we construct a schematic dissociation-metallization phase diagram and sug- 



gest experimental approaches that should significantly reduce the pressures necessary 
for the realization of the elusive metallic phase of hydrogen. 
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One of the fundamental questions relating to the properties of dense hydrogen is that 
of metallization. Since it was first predicted by Wigner and Huntington 1 in 1935, metallic 
hydrogen has long been sought after. Its value has only increased with time. The possibility 
of a high-T superconductor^, or other ordered quantum states^, has led to a tremendous 
amount of scientific investigation and thought. Despite this, many open questions remain 
as to the existence and properties of this exotic phase of matter. 

Although it was originally estimated that metallization would occur at about 25 GPa, 
pressures (P) an order of magnitude larger have since been achieved in the lab and solid 
metallic hydrogen has yet to be observed. As new methods expanded the range of pressures 
attainable by experiment, theoretical predictions for the transition grew in step. Much has 
been written about the high pressures thought to be required to obtain metallic hydrogerP^. 
In a recent experimental study, Loubeyre^ and co-workers estimate that the transition pres- 
sure (for T = K) is in the vicinity of 450 GPa. Currently these conditions are outside of 
experimental reach. 

The notion that pressures on the order of several megabar are required for metallization is 
potentially misleading, however. Indeed, metallic H has been experimentally realized already 
in the liquid phase at pressures as low as 140 GPa™ 3 . While this state was achieved for 
only a short time, it provides conclusive evidence for the existence of metallic state of H. 
Furthermore, it was recently proposed^, and subsequently confirmed experimentally^, that 
through "chemical pre-compression" hydrogen rich substances such as SiH4 (silane) metallize 
at relatively low pressures (approximately 50-60 GPa). Promisingly, silane exhibits^ a rapid 
increase in the superconducting transition temperature, T c . These works highlight the fact 
that the conventional wisdom, essentially "push harder", may not be the best approach to 
obtaining metallic H. 

In this article, we discuss the connection between molecular dissociation and metalliza- 
tion, a subject that has attracted considerable debate. Based on this consideration, we 
construct a schematic phase diagram for the liquid. We consider possible routes to the real- 
ization of the elusive phase of metallic H - ones that involve pressures lower than have been 
stated necessary, and at temperatures below what was required^ experimentally. 

While it is widely accepted that hydrogen (either as H or H2) is an insulator at low 
densities, and should be metallic at extremely high densities, it is less clear where, and by 
what mechanism, the transition between these two phases occurs. At low densities, H atoms 
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(or molecules) are sufficiently far apart that all possible atomic configurations are insulating. 
By insulating we mean that Yjp ermi falls below the Mott mobility edge^; at T = K the 
system will not conduct electricity. At extremely high densities, the opposite is true. All 
possible structural configurations are metallic - electronic states at Epermi are delocalized and 
conductivity is possible in the absence of electronic excitations. At intermediate densities, 
the situation is more complex. 
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FIG. 1. Predictions for the metallization density, r me t a i, of hydrogen based on various theoretical 
approaches. For r s > r meta i, hydrogen will be an insulator. For r s < r meta i it will be a metal. 
Predictions for bcc atomic hydrogen are shown as grey vertical dashed lines^HtLS A n estimate based 
on the Goldhammer-Herzfeld criterion for atomic (green vertical dashed line) and molecular (orange 
vertical dashed line) phases are also indicated. For comparison, a prediction for the molecular- 
atomic transition lineal based on the stability of molecular species in first principles molecular 
dynamics simulations is shown as a red solid line (dashed red line is an extrapolation). Upward 
and downward triangles bracket this transition based on molecular survival probabilities of 67% 
and 33% respectively (see^ for more discussion). For temperatures below the molecular-atomic 
transition line the liquid will be molecular, otherwise, it will be atomic. 

To frame our discussion, we begin by reviewing results for a model system: bcc atomic 
H. Depending on the approximations used, bcc atomic H (which is not thermodynamically 
stable over the conditions discussed here, i.e. at least to r s =1.30) is predicted to metallize 
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at densities as low as r s = 2.78 and as high as r s = 1.6£p^HIE] Here the density is denoted in 
terms of the Wigner-Sietz radius, defined by V/N = |7r(r s ao) 3 , where V is the volume, N 
the number of electrons, and a® is the Bohr radius. These values are plotted as grey vertical 
dashed lines in Fig. [I] If we consider values only from the more exact methods (QMC and 
GW + U), bcc atomic H is predicted to metallize within a range of densities corresponding 
to r s m 2.3 — 2.2. Achieving such a density in the liquid or solid phase requires only a 
moderate degree of pressure (P < 30 GPa, see for example^^). 

As stated above, under these conditions an atomic bcc phase is not thermodynamically 
preferred. Rather, this range of densities corresponds to phase / hydrogen - a hep lattice 
of freely rotating H2 molecules. This phase, along with phases II and III, is an insulating 
one. In fact, solid H has been shown to be an insulator at pressures up to 342 GPa by 
experiment!^ and by theory alike^. This is not a failing of the above predictions however, 
but rather illustrates a crucial point - metallization of hydrogen depends on whether the 
system is molecular or atomic. 

In general, the larger the fraction of molecules present, the higher the metallization 
density will be. For liquid hydrogen, the degree of dissociation required for metallization is 
completely determined by the density of the system. One can define a parameter r mx to 
indicate the density where a fraction X of all atoms must be unpaired in order for the system 
to metallize. For r s > r mi0 , the system will remain insulating even when all H2 molecules 
are destroyed (by T or some external means). Though subsequent heating may allow for 
conductivity by thermally excited carriers, the system will not be metallic. At most it will 
become a hot plasma. For r s < r mo , the system (liquid or solid) will always be metallic. 

The relationship between metallization and dissociation can be probed using different the- 
oretical approaches. Predictions of electrical conductivity^ 12 ^, and polarizability^ are con- 
sistent with the picture of a metallic atomic fluid for pressures corresponding to r s < r m i.o- 
An additional approach to confirming metallization is to consider the nature of interatomic 
interactions that exist within the liquid. In Fig. [2j we demonstrate the relationship be- 
tween metallization and dissociation by constructing an effective interatomic potential using 
a force-matching technique^. This approach uses configurations and forces obtained from 
first principles calculations to fit an averaged two-body potential. To generate relevant con- 
figurations, we performed large scale (as many as 1024 atoms) first principles molecular dy- 
namics simulations in the vicinity of the molecular-atomic transition line over a wide range 
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of densities. Simulations were carried out in the N, V, T ensemble, where a Nose-Hoover 
thermostat was coupled to the ionic degrees of freedom and electronic states were populated 
according to the Fermi-Dirac distribution. The electronic density (within DFT-PBE) was 
optimized at each molecular dynamics step (i.e. Born Oppenheimer MD). 
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FIG. 2. Effective two-body potentials derived (using the method oP^) from first principles molec- 
ular dynamics data. The survival probability of paired atoms, IT(t), which relates directly to the 
degree of dissociation!^, is noted next to each curve. Additionally, we note the time average value 
of the fraction of mutual nearest neighbors, < H — H >. This quantity describes the local structure 
of the liquid, but accounts poorly for changes in molecular stability. At lower densities, a), H-H 
interactions remain attractive in the dissociated (unstable) liquid, whereas at high density, b), 
they become repulsive. Small changes in the structure of the liquid (characterized by < H — H >) 
can result in significant changes in interatomic interactions and hence stability when the liquid is 
highly compressed. 

From our effective potentials, we identify a change in interatomic interactions as the liquid 
is compressed. At low densities (Fig. [2j a), H atoms are attracted to one another even when 
the system is highly dissociated (unstable). While the depth of the minima is somewhat 
reduced because of increased molecular scattering, it's overall character is essentially un- 
changed. This is consistent with an insulating atomic fluid. Conversely, for the high density 
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case (Fig. |2j b), thermal dissociation results in a marked change in interactions. The attrac- 
tive region of the potential present at low T completely disappears as the system is heated, 
even when only a vanishing fraction of molecules are destroyed (see caption of Fig. [2J. The 
resulting potential bears no resemblance to a molecular system. H-H interactions are instead 
typical of a metallic system, where ion cores interact via a screened Coulomb repulsion. This 
fact can be seen more clearly by comparing our numerical potential with one derived using 
quadratic response (method described irP*j), as in Fig. [3j This confirms not only the metallic 
character of the liquid, but suggests that such potentials, at least in the atomic regime, have 
little state dependence. We base this conclusion on the fact that our analytic potentials 
are derived without any structural information - the electronic density, which governs the 
degree of screening, is the only input variable. This implies that for the region of the phase 
diagram where the system is metallic, such potentials could reliably be used to accurately 
derive other physical properties. Furthermore, it would seem that they could be used with 
more exact methods such as Path Integral Monte Carlo, where both ions and electrons are 
treated as quantum mechanical particles. 

Based on these considerations, we can now construct a dissociation-metallization schematic 
phase diagram (Fig. [4]). In the limit of P ~ 0, molecular dissociation is governed by the 
competition of the energy cost of dissociation and the corresponding entropic gain. In this 
regime, the system undergoes a series of thermally induced transitions, from an insulating 
molecular fluid, to a semiconducting atomic fluid, until finally it becomes a conducting 
atomic plasma. In this case, conversion of H 2 — > 2H occurs continuously. 

As hydrogen undergoes compression from vanishing densities, the temperature required 
to dissociate molecules initially increases. Upon sufficient compression, however, many-body 
interactions will lower the dissociation energy and thus the temperature required for H 2 — ■> 
2H conversion. In condensed systems, the slope of the molecular-atomic transition line is 
therefore negative. 

At densities where r s < r mi this transition (H2 — > 2H) is affected by metallization. 
Metallization results in charge derealization which in turn inhibits intramolecular bonding 
through screening of the ion cores. This has the effect of destabilizing H2 bonds. For a system 
at a density of r mx , once X of the molecules have been thermally destroyed, metallization 
will occur and all remaining H 2 will dissociate. 

Since the dissociation fraction necessary for metallization decreases with density, at some 
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FIG. 3. Comparison of pair potentials generated from first principles molecular dynamics data 
(black dotted curve, r s =1.50, T=2000K) and an analytic potential derived from quadratic response 
theory (orange solid curve). The analytic potential is independent of the structure of the liquid 
and depends only on the electronic density. The region corresponding to the 1 st coordination shell 
is highlighted in purple. 

point the molecular-atomic and the insulator-metallic transition lines must merge. This 
is the point where only an infinitesimal fraction of H 2 need dissociate for the system to 
metallize. This transition should be abrupt, and may be discontinuous, i.e. 1 st order. We 
note that that the occurrence of a such a transition does not preclude the possibility of 
other structural transitions at lower pressures^. Beyond r mo o , the liquid is expected to be 
metallic for all T. A sketch indicating the location of this transition, along with all other 
phase boundaries we have discussed is given in Fig. |4| 

The fact that dissociation is not necessarily equivalent to metallization may partially 
explain the large difference in temperature between predictions of the molecular-atomic 
transition and measurements of the semiconducting-metal transition! 11 * 27 * 28 *. Indeed, Hood 
and GallP^ showed that for the case of shockwave experiments, metallization occurs when 
the fraction of H2 present is approximately 0.60, not 0.05 as was initially proposed (see 
not<P and RefP). 
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FIG. 4. Sketch of the liquid hydrogen phase diagram over a wide range of conditions. The 
density corresponding to metallization for different degrees of dissociation is indicated r mx (1.0 
and 0.5 are shown). As density increases, the degree of dissociation required for metallization 
decreases. Yellow shaded region indicates where the liquid should be metallic. Insulating-metallic 
transition is indicated by an orange line. At fixed density the liquid can sample both insulating 
and metallic configurations, therefore this transition is continuous. Black dashed line indicates 
the molecular-atomic phase boundary that would exist if metallization did not occur. T m i n is the 
minimum temperature at which the liquid is stable (either = K or the melting temperature of 
the solid phase). The intersection between the H2 — ► 2H line and the orange metallization line is 
point where the onset of dissociation and metallization would be coincident. This should result 
in a discontinuous molecular-atomic transition. We note that other first-order transitions {e.g. 
structural^) may also occur at lower densities. 

We return now to the issue of pressure and metallization. Pressurization serves several 
distinct purposes. At ambient conditions, the density of hydrogen is below the critical 
density necessary for metallization (i.e. r s > r mcta i). Application of moderate pressure 
is therefore necessary to compress the system so that r s < r meta i. Once this density has 
been achieved, metallization can be induced by two distinct mechanisms: increasing the 
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fraction of atomic H through dissociation or increasing the density of the system, thereby 
lowering the atomic fraction necessary for metallization. An increase in density also has the 
secondary effect of reducing the dissociation energy of H 2 ; covalent bonding is inhibited by 
screening that results from many-body interactions. These mechanisms may act in concert, 
as for example in a shockwave experimental. When a shockwave propagates through a 
sample, matter it encounters is compressed and heated. Compression reduces the atomic 
fraction necessary for metallization, while heating promotes dissociation due to the increase 
in molecular vibrational energy. 

Other possible routes to dissociation, and therefore metallization, include experiments 
that induce electronic or vibrational transitions (e.g. an optical probe). Work by Nagao et. 
aP^ suggests that in the solid phase, development of intermediate range crystalline order 
extends the range of pressure over which molecular phases are thermodynamically preferred. 
The corollary of this hypothesis is that at high densities, disruption of this order by way of 
chemical doping should aid in the dissociation of molecules. Similarly, metastable phases 
obtained by quenching the liquid may be more conducive to metallization. Promotion of 
collective excitations such as plasmons will also reduce the stability of molecular species. In- 
jection of hot electrons could also be used to promote molecular destabilization. Subjecting 
the system to strong electric fields or forcing charge transfer^ES may also be a viable ap- 
proaches. This could be achieved by donation from electropositive species^ or by photonic 
emission from metal surfaces with low work functions. 

By the same token, it should be possible to suppress metallization by inserting neutral 
additives (e.g. helium) into the liquid, effectively lowering the concentration of the hydrogen 
subsystem. Finally, we note that additives which catalyze H 2 — > 2H should aid in metalliza- 
tion, provided this effect is larger than the corresponding dilution. For example, transition 
metal atoms have been proposed as a means of dissociating In such cases it will be 

crucial to maintain a well defined hydrogen subsystem - metallic hydrogen is not simply 
hydrogen in the presence of a metal. 

More recently, experimental work on silane/H 2 mixtures^ indicates that under pressure, 
SiH4 has the effect of weakening H-H bonds. At the highest pressures observed, the H-H 
vibron was found to be 4090 cm -1 . As the authors point out, to soften the H-H stretch 
mode to this extent in pure H 2 requires pressures nearly 4 times as large. This result is 
encouraging, and suggests that metallic H may finally be within reach. 
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Any combination of the processes we have discussed (as well as "chemical pre-compression" ) 
should further reduce densities, and therefore pressures, required to achieve metallization. 
We therefore suggest that future high pressure experimental work focus on techniques aimed 
at promoting dissociation of H2, as it will ultimately lead to metallization. 
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